Finite-thickness effects in ground-state transitions of two-electron quantum dots 
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Using the exactly solvable excitation spectrum of two-electron quantum dots with parabolic po- 
tential, we show that the inclusion of the vertical extension of the quantum dot provides a consistent 
description of the experimental findings of Nishi et al. [Phys.Rev.B 75, 121301(R) (2007)]. We found 
that the second singlet-triplet transition in the ground state is a vanishing function of the lateral 
confinement in the three-dimensional case, while it always persists in the two-dimensional case. We 
show that a slight decrease of the lateral confinement leads to a formation of the Wigner molecule 
at low magnetic fields. 

PACS numbers: 73.21. La, 75.75.+a, 71.70.Ej, 73.23.Hk 



Two-electron quantum dots (QDs) have drawn a great 
deal of experimental and theoretical attention in recent 
years |l|, |2J, |3j- Experimental data including transport 
measurements and spin oscillations in the ground state 
under a perpendicular magnetic field in two-electron QDs 
may be explained as a result of the interplay between 
electron correlations, a two-dimensional (2D) lateral con- 
finement and magnetic field. A 2D interpretation of ex- 
periments, however, leads to inconsistencies 0, HI, pro- 
viding, for example, too low values of the magnetic field 
for the first singlet-triplet (ST) transition. There is no 
consensus on origin of this disagreement, since various 
experiments are dealing with different QDs. Evidently, 
it is important to understand basic sources of such incon- 
sistencies from view point of possible technological appli- 
cations, since QDs may provide a natural realization of 
quantum bit. This problem is also related to fundamen- 
tal aspects of strongly correlated finite systems, which are 
different from bulk and can be controlled experimentally. 

It was predicted that the ground state of N-electron 
QD becomes the spin polarized maximum density droplet 
(MDD) at high magnetic field. For a two-electron QD 
it is expected that the MDD occurs after a first ST tran- 
sition (see discussion in 0]). Theoretical calculations [f| 
assert that after the first ST transition the increase of 
the magnetic field induces several ground state transi- 
tions to higher orbital-angular and spin-angular momen- 
tum states. This issue was addressed in transport study 
of the correlated two-electron states up to 8T and 10 T 
in a lateral 0] and vertical [1, Q QDs, respectively. It is 
quite difficult to detect the structure of ground states af- 
ter the first ST transition in a lateral QD due to a strong 
suppression of the tunnel coupling between the QD and 
contacts. Altering the lateral confinement strengths, the 
transitions beyond the first ST transition are reported in 
vertical QDs [8J. In fact, the variation of the confining 
frequency with the same experimental setup opens a re- 
markable opportunity in the consistent study of effects 
of the magnetic field on electron correlations. One of 
the major aspects of the present paper is to demonstrate 
that the experimental results found in RefH can be ex- 



plained if one takes into account the 3D physical nature 
of the QD. We will discuss the additional criterion to dis- 
tinguish the 2D and 3D nature experimentally and will 
analyze the formation of the Wigner molecule in the 3D 
two-electron QD. 

Three vertical QDs with different lateral confinements 
have been studied in the experiment Q. In all sam- 
ples clear shell structure effects for an electron number 
N = 2, 6, ... at B = T have been observed, implying 
a high rotational symmetry. Although there is a suffi- 
ciently small deviation from this symmetry in sample C 
(from now on in accordance with the list of RefQ), a com- 
plete shell filling for two and six electrons was observed. 
Such a shell structure is generally associated with a 2D 
harmonic oscillator (x-y) confinement [2|. However, it is 
noteworthy that a similar shell structure is produced by 
a 3D axially-symmetric harmonic oscillator (HO) if the 
confinement in the z-direction uj z = l.buo is only slightly 
larger than the lateral confinement (lo x = w y = wq). In 
this case six electrons fill the lowest two shells with Fock- 
Darwin energy levels with n z = 0. It was found also 
that the lateral confinement frequency for the axially- 
symmetric QD decreases with the increase of the electron 
number [lfj, since the screening in the lateral plane be- 
comes stronger with large electron number. In turn, this 
effectively increases the ratio lOz/loq making the dot to be 
more " two-dimensional" , since the vertical confinement is 
fixed by the sample thickness. Indeed, the N-dependence 
of the effective lateral frequency is observed in Q. All 
these facts imply that the three-dimensional nature is a 
prerequisite of a consistent quantitative analysis of small 
QDs with a few electrons. 

Our analysis is carried out by means of the exact di- 
agonalization of the Hamiltonian for two 3D interacting 
electrons in a perpendicular magnetic field: 
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Here a = e 2 /47r£ £r and H spin — g*^s(si + s 2 )B de- 
scribes the Zeeman energy, where [Lb = eh~/2m e c is the 
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Bohr magneton. The effective mass is m* = 0.067m e , the 
relative dielectric constant of a semiconductor is e T = 12 
and \g* \ = 0.44 (bulk GaAs values). For the perpendic- 
ular magnetic field we choose the vector potential with 
gauge A = iB x r = ^B(—y,x,0). The confining po- 
tential is approximated by a 3D axially-symmetric HO 
U(r) = m*[ujQ (x 2 + y 2 ) + uj 2 z z 2 ]/2, where Huj z and Tvujq 
are the energy scales of confinement in the z-direction 
and in the rry-plane, respectively. 

The evolution of the ground-state energy of a two- 
electron QD under the perpendicular magnetic field can 
be traced by means of the additional energy Afi = 
H{2, B) - fi(l, B), where fi(N, B) = E(N, B) - E(N - 
1,B) and E(N,B) denotes the total energy of the QD 
with N electrons under a magnetic field of the strength 
B [2j. Fitting the B-field dependence of the first and 
second Coulomb oscillation peak positions to the lowest 
Fock-Darwin energy levels of the 2D HO with the poten- 
tial m*WQr 2 /2, Nishi et al. [1| estimated ujq for all three 
samples A, B, C. Although the general trend in the exper- 
imental data is well reproduced by the 2D calculations, 
the experimental positions of the ST transition points are 
systematically higher (see Fig. 3 of Ref.[l]). Different lat- 
eral confinements in the above experiment are achieved 
by the variation of the electron density, without chang- 
ing the sample thickness. Using the "experimental" val- 
ues for the lateral confinement and the confinement fre- 
quency lo z as a free parameter, we found that the value 
huj z = 8 meV provides the best fit for the positions of 
kinks in the additional energy 

Afj, = Uujqe — E(l, B) + E z (2) 

in all three samples. Here, HloqE is the relative energy, 
E(l,B) = hui + Hlu z /2 and the Zeeman energy Ez is 
zero for the singlet states (we recall that the total Hamil- 
tonian (fTJ) is separated onto the center-of-mass Hamil- 
tonian, the Hamiltonian of the relative motion and the 
Zeeman term). For the sake of illustration, we display in 
FigUthe magnetic dependence of the experimental spac- 
ing between the first and the second Coulomb oscillation 
peaks AV^, = V g (2) — V g (l) for samples A-C, which can 
be transformed to the additional energy A/z (see details 
in RefsH, 0). In the AV g — B plot, ground state transi- 
tions appear as upward kinks and shoulders @. It was 
found from the Zeeman splitting at high magnetic fields 
that \g*\ = 0.3 § and we calculate the additional energy 
with this and the bulk values. 

We nicely reproduce the experimental position of the 
first ST transitions at B — 4.2, 3, 2.3 T in samples A, B, 
and C, respectively (see Fig|T]). When the magnetic field 
is low, a difference between the calculations with different 
\g*\ factors is negligible. Upon decreasing the lateral con- 
finement Hu>o from sample A to sample C (the increase of 
the ratio lo z /u;o), the Coulomb interaction becomes dom- 
inant in the interplay between electron correlations and 
the confinement [11| ■ In turn, the smaller the lateral con- 
finement at fixed thickness is (the stronger is the electron 
correlations) the smaller the value of the magnetic field 



is at which the ST transitions or, in general, crossings 
between excited states and the ground state may occur. 
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FIG. 1: (Color online) The magnetic dependence of the ad- 
ditional energy A/i in two-electron QDs with lateral con- 
finements huo — 4.2,3.7,3.5,2.9 meV (the first, second and 
fourth values are experimental values for sample A, B, and 
C, respectively Q). The confinement in the third (z) di- 
rection hu>z = 8 meV is fixed for all samples. The results 
for \g*\ — 0.3(0.44) are connected by solid (dotted) line 
for Tiu>o — 4.2, 3.5, 2.9 meV and by dashed (dotted) line for 
hcjQ — 3.7 meV. The solid grey lines display the experimen- 
tal spacing AV g as a function of B. The arrows identify the 
position of experimental ground state transitions [||. 

There is no signature of the second crossing in the 
ground state for sample A at large B (up to 10 T). Here, 
the ratio lo z /loq rs 1.9 and the effect of the third dimen- 
sion is most visible: the confinement has a dominant role 
in the electron dynamics and very high magnetic field 
is required to observe the next transition in the ground 
state due to electron correlations. Thus, the MDD phase 
survives until very high magnetic fields {B ~ 10 T). 

A second kink is observed at B = 7 T in sample B 
Q. Our calculations with the "experimental" lateral 
confinement Tuoq — 3.7 meV produces the second kink 
at B = 9.5 T, which is located higher than the experi- 
mental value. The slight decrease of the lateral frequency 
until Tujjq = 3.5 meV shifts the second kink to B = 8.7 
T, improving the agreement with the experimental posi- 
tion of the first ST transition as well. In addition, the 
use of \g*\ = 0.3 (instead of the bulk value) with the lat- 
ter frequency creates a plateau, which bears resemblance 
to the experimental spacing AV g . However, there is no 
detailed information on this sample and we lack a full un- 
derstanding of this kink. It seems there is an additional 
mechanism responsible for the second kink in sample B. 

The most complete experimental information is related 
to sample C and we also study this sample in detail. In 
sample C the first experimental ST transition occurs at 
B = 2.3 T, while the signatures of the second and the 
third ones are observed at B « 5.8,7.1 T, respectively 
(see Figjl]). The 2D calculations (with the "experimen- 
tal" values hu> — 2.9 meV, \g*\ — 0.44) predict the first, 
second and third ST crossings at lower magnetic fields: 
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B = 1.7,4.8,5.8 T, respectively (see Fig©. The results 
can be improved to some degree with \g*\ = 0.3. To re- 
produce the data for Afj, Nishi et al. [8J have increased 
the lateral confinement (TkJo = 3.5 meV, \g*\ — 0.44). As 
a result, the first, second and third ST transitions occur 
at B = 2, 6.3, 7.5 T, respectively. Evidently, 2D calcula- 
tions overestimate the importance of the Coulomb inter- 
action. The increase of the lateral confinement weakens 
simply the electron correlations in such calculations. In 
contrast, the 3D calculations reproduce quite well the 
positions of all crossings with the " experimental" lateral 
confinement hu>o = 2.9 meV at B = 2.3,5.8,7.1 T (see 
Fig®. 
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FIG. 2: (Color online) Magnetic field dependence of the addi- 
tion energy A/i for the 2D model with hu>o = 2.9, 3.5 meV and 
for the 3D model (huio — 2.9 meV, hui z = 8 meV). Ground 
states are labeled by (Al, S), where Al and S are the total 
orbital momentum and the total spin, respectively. Grey ver- 
tical lines indicate the position of the experimental crossings 
between different ground states. 

One of the questions addressed in the experiment [8[ 
is related to a shoulderlike structure observed in a small 
range of values of the magnetic field (see our FigQ] and 
Fig. 4 of [8y). This structure is identified as the second sin- 
glet state (2, 0) that persists till the next crossing with 
the triplet state (3,1). According to Ref|8|, the ground 
state transition from the triplet (1,1) state to the singlet 
(2, 0) is associated with the collapse of MDD state for 
N = 2. Therefore, a question arises: at which conditions 
it would be possible to avoid the collapse of the MDD 
phase (in general, to preserve the spin-polarized state); 
i.e., at which conditions the singlet (2,0) state never will 
show up in the ground state. In fact, the collapse of the 
MDD depends crucially on the value of the lateral con- 
finement and the dimension of the system. We found that 
in the 2D consideration the (2, 0) state always exists for 
experimentally available lateral confinement (see FigJ3]). 
Moreover, in this range of ojq the 2D approach predicts 
the monotonic increase of the interval of values of the 
magnetic field AB, at which the second singlet state sur- 
vives, with the increase of the lateral confinement. In 
contrast, in the 3D calculations, the size of the interval is 
a vanishing function of the lateral confinement for a fixed 



thickness (Tiuj z = 8 meV). It is quite desirable, however, 
to measure this interval to draw a definite conclusion and 
we hope it will done in future. 
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FIG. 3: (Color online) The interval AB in which the singlet 
state (2, 0) survives as a function of the lateral confinement 
for 2D and 3D calculations. The confinement in the third (z) 
direction hu z — 8 meV is fixed for the 3D calculations. 

As discussed above, the decrease of the confinement 
at fixed thickness increases the dominance of the elec- 
tron correlations in the electron dynamics. Furthermore, 
this decrease, related to the decrease of the electron den- 
sity [1, [l(| , creates the favorable conditions for onset of 
electron localization. This localization (crystallization) 
in QDs is associated with the formation of the so-called 
Wigner molecule [l2j . In the 2D approach the crystalliza- 
tion is controlled by the ratio of Coulomb and confine- 
ment strengths Rw — {a/lo)/huj (l a = (Ti/m*ujQ) 1 t 2 ) 
(cf [Hj]), which is about Rw ~ 3 for the QDs consid- 
ered in experiments Q. For a 2D two-electron QD, it 
is predicted that the Wigner molecule can be formed for 
Rw ~ 200 at zero magnetic field [14|, or at very high 
magnetic field [l5| (for huo ~ 3meV and small Rw such 
as in the experiments [H). In the 3D axially-symmetric 
QDs the ratio between vertical and lateral confinements 
(anisotropy) may, however, affect the formation of the 
Wigner molecule. This problem can be analyzed by dint 
of the electron density 



i(v) = J [|*(r,r')| 2 + |*(r',r)| 2 ] dr', 



(3) 



when one electron is at the position r if the other one 
is located at a position r'. A criterion for the onset of 
the crystallization in QDs can be the appearance of a 
local electron density minimum at the center of the dot 
[l6| . For 2D QDs this leads to a radial modulation in the 
electron density, resulting in the formation of rings and 
roto- vibrational spectra |17( . 

Our analysis of the conditions realized in the experi- 
ments Q predicts very high magnetic fields (B > 12 T) 
for the formation of the Wigner molecule. However, with 
a slight decrease of the lateral confinement, at Tiluq = 2 
meV we obtain the desired result. The 3D analysis of 
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FIG. 4: (Color online) Top: the magnetic dependence of the 
ground state in 2D and 3D (hco z = 8 meV) approaches for 
a lateral confinement JkJo — 2 meV. The 2D (left) and 3D 
(right) electron densities are displayed for different ground 
states (M,S) at corresponding magnetic fields. The largest 
3D density grows from a central small core over a ring to a 
torus with the increase of the magnetic field. 



electron density (see FigQJ gives an unequivocal answer 
that at B > 7.25 T the triplet state (5, 1) can be associ- 
ated with a formation of the Wigner molecule. There is 
an evident difference between the 2D and 3D approaches: 
the 2D calculations predict the crystallization at lower 
magnetic field (AB ~ 1 T). The further increase of the 
magnetic field leads to the formation of a ring and a torus 
of maximal density in 2D- and 3D-densities, respectively. 
Notice that if geometrical differences are disregarded, 3D 
evolution of the ground state can be approximately repro- 
duced in 2D approach with the effective charge concept 
[HI (see also 0). 

Summarizing, we have shown that the confinement in 
the z direction is important ingredient for the quantita- 
tive analysis of the experimental data for two-electron 
axially-symmetric vertical QDs. In contrast to the 2D 
description, the 3D approximation provides a consistent 
description of various experimental features: the energy 
spectrum for small magnetic field, the value of the mag- 
netic field for the first and the second singlet-triplet tran- 
sitions. We propose a criterion for the additional spectra, 
which evidently demonstrates the effect of the third di- 
mension. According to this criterion the singlet state 
(2, 0) is a vanishing function of the lateral confinement 
(see FigJ3]) in the vertical magnetic field in the two- 
electron axially-symmetric vertical QD. We found that 
the decrease of the lateral confinement in the experiment 
[H until hu> = 2 meV would lead to the formation of the 
Wigner molecule at B ~ 8 T. 
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